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ABSTRACT 

We reconsider the model of quantum mechanics violation in the Kq-Kq system, 
due to Ellis, Hagelin, Nanopoulos, and Srednicki, in which CP- and CPT-violating 
signatures arise from the evolution of pure states into mixed states. We present a 
formalism for computing time-dependent asymmetries in this model and show that 
present data constrains its parameters significantly. In the future, this model will 
be put to very stringent tests at a factory. We present the theory of these tests 
and show the relation between particular decay correlations and the parameters 
of quantum mechanics violation. 
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1. Introduction 



The propagation and decays of Kq and Kq mesons has provided elementary par- 
ticle physicists with the most fruitful system for probing the fundamental discrete 
symmetries of Nature. The paradoxes of this system led Lee and Yang to postulate 
the violation of parity in the weak interactions!^' More detailed exploration led to 
the discovery of CP violation, and this system remains the only place in which 
CP violation and T violation have been observed. New experiments at CPLEAR 
are now strengthening the precision of our knowledge of neutral kaon physics, and 
beautiful new experiments on the discrete symmetries are planned for future (f) 
factories.'^^'*' 

One of the goals of current and planned experiments on the Kq-Kq system is 
to search for the violation of CPT. To interpret these experiments, one should have 
some idea of the source of CPT violation. Within the context of local quantum 
field theory, CPT conservation is a theorem!^' Thus, theories of CPT violation must 
necessarily step outside the standard assumption that particle physics is governed 
by a local quantum field theory. In the 1960's, one could plausibly build a theory 
of CPT violation in the context of a nonlocal model of the strong interactions. 
However, with the triumph of the standard model, that route is no longer open. 

On the other hand, developments in the quantum theory of gravity have opened 
another possible avenue to theories of CPT violation. Building from his results on 
the spectrum of radiation from black holes. Hawking has proposed that the gener- 
alization of quantum mechanics which encompasses gravity allows the evolution of 

ryi 

pure states into mixed states. Page then showed that any such dynamics also leads 

[^] 

to conflict with CPT conservation. These ideas raised the interesting possibility 
that one could flnd observable CPT violation due to a mechanism that lies not 
only beyond a local quantum description but also beyond quantum mechanics alto- 
gether. The notion that gravitation effects beyond quantum mechanics can affect 
elementary particle physics is controversial. For example, in theories which allow 
the evolution of pure states into mixed states, there is a serious conflict between 
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energy-momentum conservation and locality. In this paper, however, we will use a 
formalism which avoids the issue of locality. Our main goal will be to discuss how 
the issue of the violation of quantum mechanics can be addressed experimentally. 

To our knowledge, the first attempt at a phenomenological analysis of this 
type of quantum mechanics violation was made by P.H. Eberhard in the early 
70s','"' as an attempt to probe the postulates of quantum field theory. Eberhard 
suggested various tests of the existence of a unitary S-matrix. His ideas motivated 
an experimental test of quantum mechanics in the Kq-Kq system by Carithers 

[12] 

et ai. The suggestion that Hawking's idea could be tested experimentally in the 
Kq-Kq system was put forward independently by Ellis, Hagelin, Nanopoulos, and 
Srednicki (EHNS).^ 'EHNS set up an evolution equation for the Kq-Kq system in 
the space of density matrices. Their equation contains three new CPT violating 
parameters a, (3, and 7. These parameters have the dimensions of mass and might 
be expected to be of order rnj^/mpi ~ 10^^^ GeV. 

These authors were attracted to the Kq-Kq system by the fact that it con- 
tains phenomena which depend on quantum coherence over macroscopic distances. 
Though there are other phenomena that depend on macroscopic quantum corre- 
lations, for example, the persistence of superfiuid and superconducting currents, 
the Kq-Kq system gives a controlled setting, involving only one particle, in which 
precision experiments can reveal small deviations from the predictions of quantum 
mechanics. A second such system, also discussed by EHNS, is found in macro- 
scopic neutron interferometry. Here, experiments of Zeilinger, Horne, and ShuU 
have constrained a similar dimensionful parameter of quantum mechanics viola- 
tion to a level less than 0.8 x 10~^^ GeV, '"'under the assumption that quantum 
mechanics violation can randomize the neutron spin. 

Recently, Ellis, Mavromatos, and Nanopoulos (EMN) '^'^'^'*' reconsidered the 
analysis of EHNS for the Kq-Kq system and presented experimentally allowed re- 
gions for the parameters a, (3, and 7 which were consistent with nonzero values 
of the magnitude of the earher estimates. Their analysis suggested that this new 
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source of CPT violation might fully account for the observed CP violation in the 
Kq-Kq system. EMN also presented a microscopic theory of a, (3, and 7, based on 
string theory, which gives values of the size of the estimate above. While we do not 
accept their detailed microscopic arguments, we believe that their suggestion that 
quantum mechanics violation might be observable in present or planned experi- 
ments is an exciting one which deserves further consideration. This is especially 
true because the analysis of EMN was incomplete, as these authors themselves 
pointed out, in that it did not fully take into account constraints from the time 
evolution of the Kq-Kq state. 

It is the purpose of this paper to analyze the dynamics of the Kq-Kq sys- 
tem taking into account the possibility of CPT violation from mechanisms both 
within and outside quantum mechanics. We will use this framework to discuss 
the implications of past, present, and future experiments on the Kq-Kq system. 
Our formulae will include both the conventional CPT violation within quantum 
mechanics and the CPT violation outside quantum mechanics of EHNS. Our goal 
is to explain how to disentangle these two possible sources of CPT violation from 
one another, and how to distinguish both from conventional CP violation. 

Our discussion will proceed as follows: In Section 2, we will review briefly the 
conventional parametrization of CP and CPT violation within quantum mechanics 

by the parameters e and A. We will then rewrite this discussion in the language 
of density matrices. In Section 3, we will review the formalism of EHNS and 
generalize the density matrix equations to include their parameters for quantum 
mechanics violation, a, and 7. 

In Section 4, we will work out expressions for the observable quantities in 

experiments on single neutral kaons in terms of these parameters. We will show 
that the parameter a is not significantly constrained by these experiments, but 
that j3 and 7 are tightly constrained by current data. Using the classic results of 
the Carithers et a/. '^^'and CERN-Heidelberg experiments'^^ ^"'^ and recent results 
from CPLEAR, "'we determine the values (5 = (0.32 ± 0.29) x 10"^^ GeV, 7 = 
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(-0.2 ± 2.2) X 10^^^ GeV. These bounds are much stronger than those of EMN 
and lead to the simple but important conclusion that the CP violation observed 
in the Kq-Kq system is dominantly quantum mechanical in nature and of CPT 
conserving origin. The analysis of Section 4 includes provision for CPT violation 
from within quantum mechanics, but it ignores the possibility of CPT violation 
in decay matrix elements. In Section 5, we introduce this possible additional 
complication and explain how it potentially weakens our results. 

In Sections 6-8, we will discuss precision tests for CPT violation and violation 
of quantum mechanics at the future (f) factories. A (f) factory produces K meson 
pairs in an antisymmetric pure state. In ref. 4, Peccei and collaborators showed 
that, within the context of quantum mechanical CP and CPT violation, the full 
time dependence of the decays into identical states, for example, to {7r''"7r~ , 7r''"7r~}, 
is independent of the CP violation parameters and so provides an accurate mea- 
surement of the mass difference and lifetimes of the kaon eigenstates. When we 
allow for a departure from quantum mechanics, this is no longer true. We find 
an additional time-dependent oscillation which, if observed experimentally, would 
provide direct evidence of time evolution outside quantum mechanics. Section 6 
provides the general formalism for describing this effect. Section 7 gives formulae 
from which a, and 7 can be constrained independently of one another in spe- 
cific (j) factory experiments. In Section 8, we discuss various processes in which 
the two kaons decay asymmetrically. In particular, we consider the decay into 
{7r"'"7r~ , 7r°7r°}, which provides a measurement of Re e'/e, and we display correc- 
tions to the standard formulae which appear when (5 and 7 are nonzero. 
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2. Quantum mechanics of a neutral kaon beam 



In quantum mechanics, a kaon state evolves according to the action of a Hamil- 
tonian. Even if the kaon decays, its evolution is correctly described by an effective 
Hamiltonian which includes the natural width of the state. 

One conventionally writes this effective Hamiltonian as 

H = M-'-V, (2.1) 

where M and F are Hermitian 2x2 matrices acting on states in the basis {\Kq) , |^o))- 
These states are alternatively described using the basis of CP eigenstates, 

l^i) = ^{\K,) + \Kq)) ■ \K2) = -^(|i^o) - |i^o» , (2.2) 

or the basis of eigenstates of \Ks) and \Ki). 

It is easiest to express predictions for the Kq-Kq system by trading the matrix 
elements of H for parameters which express the eigenvalues and eigenvectors of this 
matrix. We will use the conventions of Maiani|^°' which express the eigenstates of 
H in terms of parameters ejvf, which is odd under CP but even under CPT, and 
A, which is odd under both CP and CPT. Explicitly, 



No 

\Ks) = ^((1 + ^S) \Ko) + (1 - es) \Ko)) 
m^^{{l + eL)\Ko)-{l-eL)\Ko)), 



(2.3) 



where 

= CM + A , ei = CM - A, (2.4) 
and Ns, are real, positive normalization factors. We write the corresponding 



6 



eigenvalues as 

ms - = (m - Am/2) - + Ar/2) 

(2.5) 

rriL - ]fL = (m + Am/2) - ^(f - Ar/2), 

so that Am and AF are positive. 

Prom this description of the kaon dynamics in terms of states, it is easy to con- 
struct the description in terms of density matrices. With the complex Hamiltonian 
(2.1), the density matrix evolves according to 

Pk{t) = e-'^-pK{Q)e^H^\ (2.6) 



The eigenmodes of this equation are 

PL = \Kl) {Kl\ 
PS = \Ks) {Ks\ 

(2.7) 

PI = \Ks) {Kl\ 
pj = \Kl) {Ks\. 

A generic kaon beam can be decomposed into these modes. If we write the expan- 
sion coefficients at time r = as ^4^, As, Aj, Aj, we find the following general 
form for the time evolution: 



PKir) = ^LPLe-r-- + Aspse-^'- + A, e'^- e"^^-- + Aj pje'^^ e+'''^- . 

(2.8) 

When px describes a pure quantum state, pi^s Pi J ^'^^ in (2-7) and the 
coefficient of the interference term Aj is correlated with A^ and As- If we write 
the pure state as pK = \K) with \K) — ql \Kl) + as \Ks), then A^ = la^p, 

— \(^s\'^ and Ai — A*j — aio^g. When pK is a mixed state, we will find that 
these properties no longer hold. 
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In the \Ki) , \K2) basis, eq. (2.2), the four eigenmodes of the density matrix 
take the form: 

(2.9) 

Here and henceforth, we normahze the eigenmodes pi so that the largest element is 
1; the normalization factors Ns, Nl in (2.3) can be absorbed into the coefficients 
Ai in (2.8). 

Any physical property of the kaon beam can be extracted from the density 
matrix by tracing with a suitable Hermitian operator. To extract the value of the 
observable V, we write 

{V)^tT[pKOr]. (2.10) 

This expression for expectation values will remain true in the generalization of 
quantum mechanics described in Section 3. In the remainder of this section, we 
will write the operators O-p associated with the most important observables of the 
neutral kaon system. 

In principle, CP and CPT violation can show up not only in the neutral kaon 
mass matrix but also in the kaon decay amplitudes. There is already some experi- 
mental evidence for a nonzero value of the parameter e' which characterizes direct 

[22 231 

CP violation in decays to two pions. ' However, the CP violation associated with 
mass mixing is much more important. In models such as the superweak model, in 
which CP violation is the result of new physics at a scale M much greater than 
m\Y, direct CP violation typically results from higher-dimension local operators 
and so is suppressed by extra powers of (mw/M). Thus, the observation of direct 
CP violation is evidence for an origin of CP violation at the weak interaction scale, 
for example, by the Kobayashi-Maskawa mechanism. 
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Although we know strictly that CPT violation cannot be the result of a local 
quantum field theory, we believe it is nevertheless reasonable to use the dimen- 
sional analysis rules of local quantum field theory to restrict possible sources of 
CPT violation. This is certainly true in a model such as that of ref. 16, in which 
violation of quantum mechanics arises from applying an unusual averaging proce- 
dure to quantum mechanical amplitudes. When we consider models of quantum 
mechanics violation, the scale M which suppresses multipoint amplitudes must 
be very large. The appearance of CPT violation within quantum mechanics will 
signal a breakdown of local field theory, and this must take place at very short 
distances compared to those probed in the LEP experiments at the The CPT 
violation outside quantum mechanics that we will describe in the next section has 
the Planck scale as its characteristic mass scale. Combining these ideas, we ex- 
pect that CPT violation is associated with processes with the minimum number 
of external particles. Thus, in our main analysis, we will include CPT violation 
only in the time development of the neutral kaon state, and we will ignore possible 
CPT- violating contributions to the decay vertices. 

Despite this argument, many phenomenological analyses include the possibility 
of CPT-violating decay amplitudes. In this case, the experimental effects of these 
CPT-violating terms cannot be unambiguously disentangled from those of quan- 
tum mechanics violation. In the standard discussion without quantum mechanics 
violation, there is a similar difficulty in disentangling CPT violation in the decay 
amplitudes from that in the K^-Ks mass matrix. In that case, one can at least 
present constraints on combinations of CPT-violating parameters. Then one can 
argue that either the individual parameters obey similar bounds or there are un- 
natural large cancellations!^"' We will present the generalization of these constraints 
to the theory with quantum mechanics violation in Section 5. 

By the same argument as that which eliminates CPT-violating decay ampli- 
tudes, we will ignore the possibihty that these new processes violate the A^* — AQ 
rule for leptonic K decays. In his review article, ref. 20, Maiani has demonstrated 
that conventional contributions to AS ^ AQ amplitudes are also negligible. Thus, 
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we will ignore 7^ AQ effects throughout this paper. 

Within this set of assumptions, we now construct the operators O-p associated 
with the semileptonic and 2-pion decays of the neutral kaons. If we ignore violations 
of the AS" = AQ rule, only the \Kq) state can decay to TT~t^v. Then, using the 
basis (l-f^i) , \K2)) of eq. (2.2), the decay rate to this final state corresponds to the 



operator 



0^+ = |a|^|Xo)(i^o| = ^( I I |, (2.11) 



Similarly, the decay rate to tt'^I v is given by 



a 



2 



0,- = \a\'\K^){K^\ = '-^\^_^ (2.12) 

In our expressions for the decay rate of neutral kaons to two pions, we will ignore 
the possibility of direct CPT violation, but we must include an allowance for direct 
CP violation. Then the amplitudes for the decay of and |^o) to tt+tt" are 

y (2.13) 

M{Ko ^ TT+TT") = Ale'^'^ + —A*2e'^\ 

v2 

We choose the Wu-Yang convention in which Aq is real. Then the decay amplitudes 
from states of definite CP are: 

M(Ki TT+TT") = V2Aoe'^° + ReAie'^' 

(2.14) 

M{K2 TT+TT") = ilmA2e'^\ 
Then, the decay rate to tt+tt~ is given, in the basis (2.2), by the operator 

/ \V2Ao + ReA2e'^\'^ (V^Aq + ReA2e-'^){ilmA2e'^) 

\ {V2A0 + ReA2e'^) (-ilmAae-*^) lilm^ae^^P 

(2.15) 

where S — 82 — Sq. One can check this result by tracing with the density matrices 
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for the eigenstates Ki and Ks, as given in eq. (2.9). We find 

tr[0+_pL] = \eL{V2AQ + ReA2e'^)+ilYiiA2e'^\^ 



(2.16) 



from which we recover, to leading order in CP violation and A7 —1/2 enhance- 
ment, the familiar result 



l^+-l = 



I — > TT+TT 



M{Ks TT+TT") 



1 



(2.17) 



The analogous argument for the neutral pion decay leads to the operator 



Coo 



|v^Ao-2ReA2e*^|2 
{V^Aq - 2ReA2e*^)(2iImA2e"*^) 



{V2Aq - 2ReA2e-*'^)(-2iImA2e*'^) 
I - 2iImA2e*'^|2 



(2.18) 

Prom this expression, in the context of purely quantum-mechanical evolution, we 
find 



I ^00 1 



M{Kl ttV 



|,^_V2^e^f. 

Aq 



(2.19) 



The quantities 77-1 — and 7700 are parametrized in terms of CP- violating observ- 
ables e and e': 



77+- = e + e 



?7oo = e - 2e' . 



(2.20) 



When we ignore all effects outside quantum mechanics and also ignore the possi- 



bility of CPT violation in decay amplitudes, we find 



[20] 



e = eL = eM - A , 



i Im742 
V5 ^0 



(2.21) 



Prom this review of the effects of purely quantum-mechanical evolution on 
the Kq-Kq system, we now procede to the case of propagation outside quantum 
mechanics. 
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3. Kaon evolution outside quantum mechanics 



We now wish to enlarge the framework of our kaon beam equations of motion 
to allow the possibility that pure states can evolve into mixed states. In this 
section, we add terms to the density matrix equation to allow this possibility. Our 
discussion will follow the arguments of EUis, Hagehn, Nanopoulos, and Srednicki 

[13] 

(EHNS). We will extend their work in providing formulae for time-dependent 
interference phenomena. 

In conventional quantum mechanics, the density matrix obeys the evolution 
equation 

t±p=[H,p]. (3.1) 
This equation guarantees the conservation of probability 

^tr[p]=0, (3.2) 

and also a set of higher identities 

d 



tr[p-] = 0, (3.3) 



which imply that the purity of the state is not changed by quantum mechanical 
evolution. 

In a two-state system such as the neutral kaon system, the density matrix can 
be expanded 

p = pOl + pV, (3.4) 

where i — 0,1,2,3 and cr* is a Pauli sigma matrix. The statement of conservation 
of probability tr[p] = 1 becomes — ^, and the statement that p has positive 
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eigenvalues becomes 

ip'f > hp')" (3-5) 
1=1 

If the Hamiltonian is expanded in the same way, the equation of motion can be 
written 

= 2i^^WfPa^. (3.6) 

or 

The quantum mechanical description of the neutral kaon system is only slightly 
more complicated. Here we work with the non-Hermitian Hamiltonian (2.1) which 
includes a provision for the kaon states to decay. The evolution equation (2.6) 
leads to the equation 

-^p = 2e'^^Mya^ - p - r(pV^ + p'l). (3.7) 
otr 

The value of p° = 2 tr[p] will now decrease. However, the inequality (3.5) must still 
be satisfied, and it will be as long as the matrix F has two positive eigenvalues. 

Hawking proposed that the modifications of quantum-mechanical evolution 
due to quantum gravity effects could be described by writing a more general linear 
equation for the density matrix. The most general term that we could add to (3.7) 
is 

-/i°Vl - h^^'y^ - h'^'^V (3.8) 

There are two natural restrictions on these terms. First, they should be consistent 
with probability conservation. Second, they should not decrease the entropy of 
the density matrix; though pure states can evolve into mixed states, mixed states 
should not evolve into pure states. The first of these requirements sets h^^ = 0. 
The second requirement eliminates any h^'^ terms, since these would order the 
completely random distribution with p — hi, and also requires that the remaining 
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submatrix h^^ be positive definite!^'^' This leads to the following set of equations 
for the components of the density matrix: 

^; (3.9) 

dr 

where, from here on, i,j, k = 1,2,3 only. Notice that the antisymmetric part of h''^ 
can be absorbed into , so that we may assume from here on that h is symmetric. 

EHNS simplify this formalism by imposing one further assumption. If the 
action of quantum gravity is universal among flavors, the new term cannot change 
strangeness. Alternatively, if the basis chosen by the quantum gravity interactions 
is close to the basis of quark mass eigenstates, the strangeness nonconservation 
will be proportional to the square of the rotation angle between these two bases. 
If this angle is of order the Cabibbo angle, strangeness violation will be a higher- 
order effect. Under cither assumption, we may concentrate on the part of h which 
conserves strangeness. Since strangeness is measured by the operator Os = —cr^ 
in the basis of CP eigenstates, the restricted SH must satisfy 

tr[(7^/i*VV] = , i.e., h^^ = 0. (3.10) 

This leaves 

/O 0\ 

(3.11) 

7/ 

where the EHNS parameters a, P, j satisfy 

a,7>0, a7 > /3l (3.12) 

The equations (3.9) for the density matrix components are linear equations 
which can be solved by diagonahzing a 4 x 4 matrix. Since CP violation is a small 













/i = 2 





a 






U 




7 
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effect, of order 10~^ of Am and AF, a perturbative solution is quite appropriate. 
To obtain tliis solution, we first rewrite the equations (3.9) in the basis of matrix 
elements of the density matrix: 



P 




(3.13) 



It is useful to first write the purely quantum- mechanical equation (3.7) in this 
basis: 



d 



fpi\ 

P2 
P 

\PJ 



-r + 



/-Ar/2 +ield* 

Ar/2 +iesd 

-i€*gd* -ieid +iAm 

\ +iesd +ield* 



-ie^d \ 
-ie*gd* 


-iAm ) 



P 

\P ) 



, (3.14) 



with 



d = Am + -Ar. 



(3.15) 



It is straightforward to check that the parameters are chosen so that eigenvectors of 
the matrix reproduce (2.9), with the correct eigenvalues. The perturbation (3.11) 
adds to the quantity in brackets the matrix 



/-7 


7 






7 


-7 




-i(3 






—a 


a 


V-i/3 




a 


—a ) 



(3.16) 



Treating this addition to the equations in perturbation theory, and working to 
the order of the leading contribution to each matrix element, we can work out the 
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new eigenmodes and eigenvalues. The eigenmodes are: 



PL = 



PS = 



PI = 



PI 



|ez,|2 + 7/Ar + 4/3(Am/Ar)Im[eL/(i*] - /^VMP cl + (^/d 
el + P/d* 1 

1 e| - p/d* 

es - (3/d \es\^ - 7/Ar - 4(3{Am/Ar)lm[es/d*] - (3^/\d\^ 

el - (3/d* 1 
—ia/2Am es + P/d 

eL — (5/d ia/2Am 
1 + (3/d* 



(3.17) 



The corresponding eigenvalues are corrected by the shifts 
f^f + a, Am^ Am-(l-i(/3/Ar)2) 



(3.18) 



The shifts of T^,, Vs, and Am can be absorbed by redefinition of these parameters. 
The shift of Am is of relative size 10~^ and so is negligible in any event. If we 
redefine f to be the average of the new Ts and T^, then the interference terms pi 
and pj fall off at the rate 

r + (a - 7) . (3.19) 

This correction is not relevant to current experiments unless a is as large as lO^^F; 
in that case a would be 10 times larger than the familiar CP- violating parameters. 
We will retain this shift, for completeness, in our formulae below, but we will ignore 
it in our analysis of the present experimental situation. 

To summarize, under the infiuence of quantum mechanics violation as described 
by the formalism of EHNS, the most general initial density matrix evolves according 
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to 

PK{r) - ^LPLe"^^^ + Aspse-^'^ 

+ Ai PI e-(r+«-7)T g-iAmr ^ p_ g-(f +a-7)r g+iAmr ^ 

(3.20) 

where now the eigenmodes are given by (3.17). 

To understand the changes that have been induced in the evolution of density 
matrices, it is useful to rewrite the eigenmodes slightly in order to emphasize their 
similarity to (2.9). It is very convenient to introduce the split e parameters: 

ef = ei ± P/d , e| = 65 ± P/d . (3.21) 

Then the expressions for the eigenmodes can be rewritten as follows: 



PL 



PS 




+ 7/Ar + A{(3/i^T)lm[eld/d*] e+ 



4* 



es |e+|2 - 7/Ar - 4(/3/Ar)Im[e+d/d* 



(3.22) 



C J 



PI = 



Pi 



—ia/2Am 

ia/2Am 



The eigenmodes are written here to sufficient accuracy for the analysis in the rest 
of this paper. The complete expressions for these eigenmodes to second order in 
small parameters are given in the Appendix. 

Except for their smallest matrix elements, of which only pn and ps2 arc im- 
portant in practice, the expressions (3.22) take precisely the form of the density 
matrices of pure states given in (2.9). However, the various modes contain differ- 
ently shifted e parameters. We will see in the next section that different neutral 
kaon experiments are sensitive to different choices of these parameters, in such a 
way that the various sources of CPT violation can be distinguished. 
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4. Experimental determination of f3 and 7 



In this section, we will show how present experimental data constrains the 
EHNS parameters. We will show how to combine the very accurate measurements 
of the parameters from the experiments of the early 1970's with new data on 
time-dependent evolution from the CPLEAR experiment. This comparison will 
give stringent constraints on (3 and 7 which limit their effects to be at most about 
10% of the observed CP violation in the neutral kaon system. This rules out the 
possibility, suggested in ref. 16, that a deviation from quantum mechanics of the 
nature considered here is the major source of observed CP violation. We will also 
recommend a parametrization of the future, more accurate CPLEAR data which 
will facilitate comparison of this data with the EHNS formalism. 

In this section, we will consider specifically the following two observables of 
the neutral kaon system: 



5{r) 



N{K{t = 0) ^ TT+TT-) 

N{K{t) TT-i+u) - N{K{t) n+e-u) 
N{K{t) TT-i+u) + N{K{t) 7r+£-i?) 



(4.1) 



These quantities are given in terms of density matrices and the decay operators 
defined in Section 2 by 



In our analysis, we may ignore the effects of Imyl2 in 0-\ , since e'/e ~ lO"'^. 

Thus, the formulae we present for the charged pion decay apply equally well to the 
neutral pion decay. It is easy to restore the effect of e' by evaluating (4.2) more 
exactly. 
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In a K'^ beam which has evolved to large r, the quantities 5 and R-\ reflect 

the pure K]^ eigenstate. In this case, we can evaluate (4.2) with the density matrix 
PL and obtain the following results: 

Sl = 2Ree+ 

(4.3) 

Rl = PLi = + 7/Ar + 4{(3/Ar)lm[eld/d*] 

Alternatively, we could consider the evolution of a pure or X° state created 
in a strong interaction process. For state which is initially pure K^, the evolving 
density matrix is given by (3.20) with — As — Aj — Aj = ^, up to corrections 
of order e, a. (For an initial reverse the sign of the interference terms.) Then 
the time-dependent quantities (4.1) are given by* 



2 cos(AmT)e-(r+"-^)^ + 2Ree^e-^«^ + 2Ree+e-^^^ 



(4.4) 



R+-{t) = e"^^^ + i?Le-^^^ + 2|?7+_| cos(AmT + 0+_)e-(^+'^-T)^ 
with Rl as in (4.3), and 

\fl+-\e't>+- ^ el . (4.5) 

Notice that the interference term in the time-dependent asymmetry depends 
on e^, while the Kl lepton asymmetry depends on ej. Thus, by comparing time- 
dependent and long-time measurements, we can find a constraint on /5 which is 
independent of other sources of CPT violation. By comparing the value of |e^P 
to a determination of Rl, and using this bound on /3, we can also obtain a bound 
on 7. 

The difficulty in implementing this program is that experiments on the time- 
dependent evolution of neutral kaon beams typically report fits of their data to 
the conventional CP T-conserving theory, in which e^, ej, and the square root of 



★ This formula includes all corrections linear in a. 
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pii arc not distinguished. Thus, it is important to understand which particular 
parameters are actually constrained by each given experiment. 

The value of Sl poses no difficulty. The current world average'^' gives 

5l = 2ReeJ = (3.27 ±0.12) X 10"^ . (4.6) 

In the early 1970's, the CERN-Hcidclberg collaboration carried out beautiful 
studies of the time-dependence of K^-Ks interference. These studies con- 
firmed that the conventional, CPT-conserving parametrization of the neutral kaon 
system indeed gives a good description of its time-dependent phenomena. How- 
ever, it is very difficult to reconstruct the constraints that these experiments put 
on more general models of time-dependence. The experiments involved kaon pro- 
duction from a platinum target; in this situation, the proportion of versus 

initial states is not known a priori and is also momentum- dependent. The 
particle-antiparticle asymmetry must be described by an unknown function 

N{p) + N{p) ^ ' 

In ref. 18, data on neutral kaon decays to tt'^tt~ was binned in momentum and 

decay time and these distributions were used to fit for the parameters |r]_| |, 0_| , 

and A{p) (in each bin) from the assumed relation 

i?^_(r) = e"^^^ + |?7+_|2e-^^^ + 2A{p)\r]+-\ cos(Amr + 0+_)e~^^. (4.8) 

Unfortunately, the absolute magnitude of the interference term can be absorbed 
into the parameters A{p), so it is not possible to determine the coefficients Rl 

and \fj-\ I separately. Since the strongest constraint on \ri-\ | from this data set 

comes from the long-time tail of the decay distribution, we will consider the CERN- 
Heidelberg determination of |r;_| | to be a measurement of R^: 



'Rl = (2.30 ± 0.035) x 10"^ (4.9) 
Prom Fig. 4 of ref. 18, it is clear that quantum mechanics does correctly describe 
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the interference region to few-percent accuracy after adjustment of the A{p), but 
we will not make use that information in the discussion below. 

Shortly afterward, the experimental group of Carithers et al, working at the 
LBL Bevatron, used their data on X — > tttt decays in a regenerated beam 
to directly measure the magnitude of the interference term in neutral kaon evolu- 
tion. In the notation of eq. (4.4), they obtained the result 

= 0.972 ±0.021 

Very recently, the CPLEAR collaboration has pubhshed as its first result a de- 
termination of \fj-\ — I using hadronically produced neutral kaons tagged by an ac- 

[3] . 

companying charged kaon. Their determination is dominated by the interference 
region, and so we may interpret their result as a measurement of |f/_| |. This gives 

|77+_| = (2.32 ± 0.14) X 10-3 . (4.11) 

Averaging these measurements using the CERN-Heidelberg value for Rl, we find 

|77+_ I = (2.249 ± 0.054) x 10-^ . (4. 12) 



(4.10) 



Finally, we see no difficulty in accepting the world average of measurements of 
the interference phase 0-| as a determination of the phase of : 

0+_ = (46.5 ± 1.2)° (4.13) 



The comparison of these four numbers allows us to constrain (3 and 7. We can 
find /3 from the relation 

d 



Re^ = Ree;^ - Ree;^ . (4.14) 
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To analyze this relation, we use the fact that (3 is real and d is very well known: 

d = Am + ^Ar = ((3.522 ± 0.016) + i(3.682 ± 0.008)) x 10"^^GeV (4.15) 
It is convenient to parametrize 



where 05^1^ = (43.73±0.15)° is the sii^erwea/s^j/iaseof Maiani's conventions!^"' With 



this definition 

2 sin (\)sw 2 



/3 - TT-^d - cos0+_) (4.17) 



This yields 

/? = (3.2 ± 2.9) X 10-^^ GeV . (4.18) 
With this determination of /3, we can evaluate 7 by comparing Ri^ to e^: 

7 = Ar[i?i - - 4(/3/Ar)|77+_| sin(205W/ - 0+_)] . (4.19) 

which gives 

7 = (-0.2 ± 2.2) X 10"^^ GeV . (4.20) 

These constraints are similar to those obtained by Ellis, Mavromatos, and Nanopou- 
^^^[15,16] one-parameter fits for /3 and 7. However, our constraints hold in the 
general three-parameter EHNS model; they exclude the possibility that /3 and 7 
could give large individual contributions which cancel in Rj^. 

The geometry of these constraints is shown in Fig. 1. In Fig. 1(a), we show the 
systematics of the various values of e. This figure should be compared to Fig. 1 of 
ref. 15. The parameter e£, considered as a vector in the complex plane, is directly 
determined by the interference measurement. Prom the endpoint of this vector, 
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is found by making an excursion downward at 45° by a displacement proportional 
to Similarly, Ri is found by moving outward a distance proportional to 7, after 
correcting for a (3 effect. In Fig. 1(b), we show the constraints on these excursions 
given by the experimental values (4.6), (4.9), (4.12), (4.13). The various constraints 
overlap in the e plane in such a way as to constraint (3 and 7 to contribute only a 
very small part of the CP violation phenomenology. 

It is important to note that the constraints on CPT violation outside of quan- 
tum mechanics that we have considered here are quite independent of the possibility 
of CPT violation within quantum mechanics. Such a source of CPT violation can 
be constrained, just as in the analysis without quantum mechanics violation,'^"' by 
verifying the extend to which ei^ — + P/d is parallel to d in the complex plane. 
We will give the precise argument in the next section. Since we have not used 
any information on the absolute phase of e^;, in order to bound (3 and 7, the stan- 
dard constraints on CPT violation within quantum mechanics are not significantly 
weakened when we allow for the presence of the EHNS parameters. In the fu- 
ture, the measurement of will allow a stronger constraint on this type of CPT 
violation. 

To derive further constraints on the violation of quantum mechanics, we strong- 
ly recommend that time dependent decay distributions be fit directly to the for- 
mulae (4.4), to determine the five independent parameters Rl, Sl, I77-1 — |, 0-| — , a. 
This will make it possible to check, without the ambiguities of our analysis here, 
whether the rapport among these parameters definitively excludes the presence of 
new terms in the density matrix evolution equation. 



23 



5. Effects of CPT violation in decay amplitudes 



At the end of Section 2, we argued that, whether CPT violation arises from 
the breakdown of quantum mechanics or from the breakdown of local quantum 
field theory within quantum mechanics, the effects of direct CPT violation in de- 
cay ampUtudes should be small compared to those in neutral kaon propagation. 
Nevertheless, much emphasis is given in the htcraturc on CPT violation to disen- 
tangling effects of CPT violation in decay vertices from CPT violation in the kaon 
mass matrix. In this section, we will show that this separation can still be made, 
and constraints on (3 and 7 deduced, if one allows for CPT- violating additions to 
the kaon decay vertices. However, this analysis will require additional assumptions 
which, though not unreasonable, are not airtight. 

To begin the analysis, we need the more comphcated forms of the decay opera- 
tors O-p which allow for CPT violation in the decay amplitudes. The paramctriza- 
tion given in Maiani's review article'^"' leads to the following expressions, which 
replace (2.11), (2.12), (2.15), and (2.18): for the leptonic decay amplitudes. 



6p 



for the TTTT decay amplitudes, 

1 




6p 



(5.1) 




Coo = l-'^ool 




X 



Y 





K2) 


(tttt 


Ki) 



(5.2) 



(5.3) 



More explicitly. 



y+- 

>00 



(^) - 2e' 



(5.4) 



The quantities Re(So/Ao) and Re(6/a) parametrize CPT- violating decay ampli- 



24 



tudes; e' has the CPT conserving value (2.21) shifted by the amount e*'^(Rei?o/^0~ 
Rei?2/^2) / \/2 which accounts for a different degree of CPT violation in the isospin 
7 = and 1 — 2 pion-decay channels!^"' With this generalization, the value of Sl 



becomes 



5L = 2Re[e+ + ^] (5.5) 



and the parameters Rl, \rj-\ — |, and 4>-\ — in (4.4) are shifted to 

Rl = 7/Ar + |f^+_|2 + 4(/3/Ar)Im[?7+_d/d* - Y+_] 

|r7+_|e*^+-=67+F+_. 



(5.6) 



When we include these corrections into the relation between ej^ and ej, we 
find, instead of (4.17), the following expression for /3: 

pH Re 1-) = — 77+- cos0+_) . (5.7) 

^ 2sin05VF Aq' 2sin05w.^2 ' ' ^ ' 

The relation for 7 remains 

7 = ^T[Rl - |r;+-|2 - 4(/5/Ar)|r;+_| sin(205W^ - 0+_)] (5.8) 

if we ignore a very small correction proportional to (5 -IvnY^ = /3 • Ime'. Thus, 

our previous constraints on /3 and 7 now appear as constraints on combinations of 
CPT-violating parameters. 

For completeness, we should add one further constraint on a combination of 
parameters characterizing CPT violation within and outside of quantum mechan- 
ics. This constraint, which is reviewed in ref. 20, uses unitarity and the dominance 
of the isospin tttt decay channel to determine the phases of the CPT-conserving 
and CPT-violating mixing parameters. If we denote these components of the kaon 
mixing parameters as em and A, as in eq. (2.4), then this constraint implies that 
the phase of eu is (psw of ^q. (4.16), and the phase of the combination 

A -5^ (5.9) 

is (t)sw ± 7r/2. 
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To apply this relation, rewrite the second of eqs. (5.6) as: 



|^,_|e*-=.M-(A-^)-^. (5.10) 

Notice that, because d points at 45°, the last term has a phase of (—45°), so that 

it is also perpendicular to em- On the other hand, the experimental value of 0-| 

is very close to that of (t)sw-i though there is a small discrepancy: 



(/.+_ - <i)sw = (2.8 ± 1.2)° . (5.11) 

Thus, by comparing the components of the right and left hand sides of (5.10) per- 
perdicular to em, we obtain a third constraint on the parameters of CPT violation. 

This additional constraint reads: 

RcBq 



(5±\d\ 



(-5.6 ± 2.5) X 10"^^ GeV . (5.12) 



This must be combined with the results of eqs. (5.7) and (5.8): 



P + ^Re(- - ^) = (3.2 ± 2.9) x IQ-^^ GeV 

7 - 2|77+_d|Re(- - ^) = (-0.2 ± 2.2) x lO-^i GeV . 

a Aq' 



(5.13) 



In the last relation, we have used the approximation ^sw ~ 0+-- These three 
equations provide three constraints on four parameters and so cannot disprove the 

existence of CPT violation. However, they imply that, unless there arc unnatural 
cancellations among these parameters, the magnitude of CPT violation should be 
at most about a tenth that of CP violation. 
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6. Tests of quantum mechanics at a factory: formalism 



A high-luminosity factory has been recognized as a facility which gives par- 
ticularly incisive tests of CPT violation. In models in which CPT violation arises 
within quantum mechanics, Peccei and collaborators have shown how, by study- 
ing the full set of possible time- dependent asymmetries observable at a factory, 
one may disentangle CPT-violating terms in the neutral kaon mass matrix from 
CPT-violating contributions to kaon decay amplitudes!'''^^' This analysis is made 
possible by the very simple time-dependence predicted for the Kq-Kq state which 
evolves from the decay of the 0. Since the has spin 1, its decay to two spinless 
bosons produces an antisymmetric spatial wavefunction. This means that, when 
the (f) decays to two neutral kaons, those particles must remain in opposite mass 
eigenstates until one decays. This strong constraint from quantum-mechanical 
coherence governs the whole phenomenology of factory measurements. 

Because of the importance of quantum-mechanical coherence in decays, a 
factory is also an ideal place to search for terms in the kaon evolution which violate 
quantum coherence. In this section, we will explain how the standard formulae for 
kaon correlations in decay are changed by the introduction of quantum mechanics 
violation according to the model of EHNS, and we will point out particularly 
incisive measurements for determining or constraining the EHNS parameters. 

The key to our analysis will be the form of the density matrix for the two- 
particle system which results from decay. We will first review the form of this 
density matrix in the case in which CPT is violated only by corrections within 
quantum mechanics. Then we will show how this result is generalized in the EHNS 
model. 

In the following discussion, we will assume that the resonance is pure spin 1, 
with no quantum mechanics violation in its decay amplitudes. Throughout our 
analysis, we will ignore background processes which produce KqKq in a spin 
combination, and also effects of finite detector size. The influence of these effects 
in conventional factory analyses are described in ref. 5. 
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A spin 1 (j) meson decays to an antisymmetric state of two kaons. If the 
kaons are neutral, we can describe the resulting wavefunction, in the basis of CP 
eigenstates \Ki), 1-^2), as 

(f)^ -^(\Ki,p> ®\K2,-p> -\K2,p> ®\Ki,-p>^ . (6.1) 

The two-kaon density matrix resulting from this decay is a 4 x 4 matrix. We can 

express this matrix concisely by introducing the set of 2 x 2 matrices appropriate 
to the \Ki), \K2) basis: 







Then the state (6.1) corresponds to the density matrix 

P ^^(^pi® p2 + P2® Pi- P+® P- - P-® P+) ■ (6-3) 

To work with the density matrix (6.3), we must express its components pi 
in terms of eigenstates of the single-particle density matrix evolution equation. 
Then we can assign each eigenstate its natural time-dependence. This procedure 
accounts correctly the full time-dependence of correlations. 

To see how this works, we will first express P in terms of eigenmodes in the 
quantum-mechanical case. It is straightforward to express the components (6.2) in 
terms of the four matrices (2.9) and then to express the combination (6.3) in terms 
of these elements. One finds, to no great surprise, an antisymmetric combination 
of \Ki,) and \Ks). Here and for the rest of this section, we will denote results 
derived assuming quantum mechanical evolution (but not CPT symmetry) with a 
superscript diamond. Thus, 

^(^) = 1 + 2Re{eseL) ^ ^ _ ^ p. _ p. ^ p^^ _ (6 4) 

The prefactor, which we have written to order e^, corrects for the fact that \Ki) 
and \Ks) are not orthogonal. Supplying the proper time-dependence, we find the 
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density matrix for processes in which the first kaon decays at proper time ri and 
the second at proper time T2. We find 

-PI® p^-e-»^^(^i-^2)g-r(ri+r2) ^^^+iAm(n-r2)g-r(ri+r2)^ _ 

(6.5) 

It is equally straightforward to perform this computation when quantum me- 
chanics violation is included. One must first work out expressions for the compo- 
nents (6.2) in terms of the eigenmodes (3.22). To first order in small quantities, 
we find 

Pi = PS + {ll^'^)pL - eg* PI - eg Pi 
P2 = -(7/Ar)p5 + PL - e^pi - el* Pi 
P+ = -€l*PS - ^sPL + PI + i{a/2Am)pi 
P- = -(^lPs - (^s*PL - iia/2Am)pi + pj . 
Inserting these expressions into (6.4), we find 



(6.6) 



2 



(6.7) 



Ps^PL + PL® PS - PI ® Pi - Pi® PI 

- 2^(P5 ®pi + PI® PS) - ^-^\Ps ®Pi + Pi® PS) 

+ ®Pi + Pi® PL) + ^-^{PL ®Pi + Pi® PL) 

. a , X 27 . X 

+ ^~^\Pi ®Pi- Pi ® Pi) + -^[PL® PL- PS ® PS) 

In eq. (6.7), each coefficient is given to first order in small quantities. We will carry 
out our analysis to this order. However, if some of the EHNS parameters are more 

severely constrained than others, one might need the expressions of second order 
in the less constrained parameters to determine the correct limits for the more 
constrained ones. For example, we saw in Section 4 that the terms quadratic in /3 
affect the limits on 7. In case a similar situation arises in factory experiments, 
we have given in the Appendix the complete formula for P correct to second order 
in small quantities. 
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When we study time-dependent correlations, each term of (6.7) will lead to 
a characteristic exponential behavior. The first line of (6.7) is identical (within 
the approximation given) to (6.4) and thus contains only time-dependences al- 
lowed within quantum mechanics. The next three lines of (6.7), however, contain 
completely new structures. The second line of (6.7) leads to terms in the decay 
distribution which behave as 

cos(AmTi - </,)e-(^+«-7)ng-rsr2 ^ ^g g^ 

and, in the same way, with 1 and 2 interchanged. The third line leads to similar 
expressions with the decay rate F^. The first term in the fourth line leads to the 
even more bizarre time dependence 

sin(Am(Ti + T2) - 0')e-(r+«-7)(ri+r2) _ (g_9) 

If we interpret frequencies as energies, both of the forms (6.8) and (6.9) signal 
a switch from positive to negative values of the energy. This is a rather subtle 
breakdown of energy conservation, which should be expected in the framework of 
density matrix evolution equations according to the analysis of ref. 10. This subtle 
effect does not obviously lead to macroscopic violations of energy conservation. 
However, in the (p factory experiments, one does not need to wait for the problems 
of energy conservation to built up to a macroscopic violation; one can instead track 
these violations directly in the frequency dependence of corrections. Finally, the 
second term in the fourth line contains the time-dependences 

g-rs(ri+r2) ^ ^-ri(n+r2) _ ^g_;LO) 

Both terms in this line signal a breakdown of the antisymmetry of the final state 
wave function, which corresponds to a subtle breakdown of angular momentum 
conservation. 
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The basic observables computed from P are double differential decay rates, 
the probabilities that the kaon with momentum p decays into the final state /i at 
proper time ri while the kaon with momentum (— p) decays to the final state /2 
at proper time T2. We denote this quantity as n; /2, T2). If we denote the 

expression (6.7) schematically as 

P = ^AyPi®p^- , (6.11) 

where 2, j run over S,L,I,I, and write the corresponding eigenvalues as Aj, then 
the double decay rate is given by 

nfh n; /2, T2) = X; tr[piOf,] tT[pjOf,]e-^^^'-^^'-^ . (6.12) 

Since it is easier to understand a distribution in one variable, much of the 
analysis of factory experiments has made use of the integrated distribution at 
fixed time interval At — ti — T2- We will assume, in working with this quantity, 
that At > 0. Then this time interval distribution is defined as 

00 

P(/i;/2;Ar) = J d{n + T2)V{h,n; f2,T2) . (6.13) 

At 

This time interval distribution is very useful for obtaining the standard CP viola- 
tion parameters of the neutral kaon system. However, since (6.13) integrates out 
one of the exponentials in (6.12), this integral does not possess the strange time de- 
pendences which signal quantum mechanics violation. A different quantity which 
shows these unusual effects more clearly is the double decay rate interpolated to 
equal times: 

Q(/i;/2;t)= Wi,t;/2,t) . (6.14) 

In the case in which quantum mechanics is exact and the density matrix is given 
by (6.5), this expression is a linear combination of decreasing exponentials, with no 



31 



oscillatory terms. However, when we introduce the EHNS parameters, Q{r) can 
acquire terms with the oscillatory dependences cos(AmAT — 0) and cos(2AmAT — 

00 ■ 

The most striking example of this modification of the quantum mechanical 
prediction appears in the case of decay to identical final states fi — f2 — f- In this 
case, the quantum mechanical prediction for the double decay rate is especially 

simple!^' Using 

tr[psOf] tr[pLOf] = tr[piOf] tr[pjOf] = \{f \ Ks)\\f \ Kl)\' , (6.15) 
we can rewrite the expression for the double decay rate as 

V^^\f,ri;f,r2) = Cx [e-^^^i-^^^^ + e-^^^i-^«^^-2cos(Am(Ti-T2))e-^(^^+^^)] . 

(6.16) 

This quantity depends on the two times in a manner completely fixed by quantum 
mechanics irrespective of the properties of the decay amplitudes. In particular, at 
Tl = T2, the double distribution vanishes, as a consequence of the antisymmetry 
of the final state wavefunction. All of these conclusions hold whether or not CPT 
symmetry is preserved. 

On the other hand, when quantum mechanics is violated, decays to identical 
final states can have a much less constrained structure which includes all of the 
time dependences found in the most general case. We will see examples of this in 
the next section. 
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7. Tests of quantum mechanics at a factory: identical final states 



Now that we have clarified the general form of the effects of quantum mechan- 
ics violation in the EHNS model, we will present expressions for the dependence 
of particular observables on the EHNS parameters. In particular, we will show 
explicitly how the terms of the equal time distributions Q(/;/;t), defined in eq. 
(6.14), can be used to constrain these parameters. In this section, we will return 
to the framework of Section 4, in which CPT violation in decay vertices and e' are 
neglected. However, in contrast to the case of experiments on single kaons, the con- 
straints on the EHNS parameters from factory experiments are not essentially 
affected by the inclusion of CPT violation in decay vertices. In the next section, 
we will explain this point and also discuss some aspects of the measurement of 
these decay parameters. 

Consider first the case in which both kaons decay semileptonically. The cases 
in which both kaons decay to T:~t^u or to 7r"'"£~P are examples of decays to iden- 
tical final states whose special properties were discussed at the end of the previous 
section. It is straightforward to work out the double time distribution for these 
cases by using the expression (6.7), the explicit forms of the density matrix eigen- 
modes (3.22), and the decay operators (2.11), (2.12). We find, to first order in 
small parameters 



P(£±,ri;£±,T2) 



a 



X I (1 ± 4Re em) + e'^^^^-^'^^ - 2 cos(Am(ri - r2)e-(^+°-^)(^i+^2)] 

± 4-^sin(AmTi - (t)sw)e~^^^'^~^^^'e-^'^^ + (1^2) 

\a\ 

± 4^sin(AmTi + (t>swy~^^^''~^^'''e-^^^^ + (1^2) 
\d\ 

+ 2-^ sin Am(ri + rs) e-(r+"-^)(^i+^2) + 2-^ Te-r^n+r^) _ g-rs(n+r2)i \ 
Am ^ ^ Ar^ J J 

(7.1) 

Notice that the first term in the brackets has a form quite close to the canoni- 
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cal form (6.16) predicted by quantum mechanics, while the remaining terms give 
systematic corrections to this result. For comparison, 



P(£±,Ti;£T,T2) = 



Ia|4 



X I [e-^sn-Ti^r^ ^ ^_r^^,_rsr2 ^ 2 cos(Am(ri - T2))e-'^^'''+^'^] (7.2) 
+ 0(e,a,/5,7)j. 



The complete expression is given in Appendix B. 

The form of the corrections to quantum mechanics are easiest to see by inter- 
polating to the line ti — T2. On this hue, (7.2) becomes 



P(£±,T;£T,r) = ^-46-2^^ , (7.3) 



8 

and (7.1) has a similar, though less dramatic, simplification. Then one finds 

l[l-e-2("--)^(l-^sin2Amr)] 

+lJLre+Arr_ -Arr] (7-4) 
^2ArL^ J 

±2^ rsin(AmT - 05W^)e-^r^/2 ^ sin(AmT + (l)sw)e+^^^^^] 

\d\ 

The three coefficients a, /3, and 7 are selected by terms which are monotonic 
in T, oscillatory with frequency Am, and oscillatory with frequency 2 Am. 

It is amusing to note that the positivity of the expression (7.4) under the 
conditions (3.12) is maintained by a delicate interplay of the correction terms. 
This is most easily seem by examining the limits a ^ 7 > and 7 ^ a > 0. 

Similar information is provided by the decay distribution to tttt final states on 
both sides of the (f) decay process. Since we ignore e' effects, the decay distributions 
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to TT^TT and final states are identical up to the overall normalization. For 

the specific case of 7t~^7t~ decays on both sides, we find, to second order in small 
parameters 

P(7r+7r-,Ti;7r+7r-,r2) = 2|Ao|'^ 

- 2|?7+_|2cos(Am(Ti - r2))e-(^+"-^)(^i+^2) 



X 



+ 4:-^\v+-\ sin(AmTi + 0+_ - (l)sw)e~^^^''~^^^'e-^'^' + (1^2) 

\a\ 

-rs(Tl+T2) 



e 



Ar \d\ COS05VK 

(7.5) 

where |f/_| \e'^'t'+- = e~ and Rl is defined as in (4.3). Specializing to the line 

Ti = T2 = T, we find 



-ATt , 



X 



f|-^_|2N _g-2(a-7)rj 2 ^ |r/^_| ~ ^'^'^^^ ^ ^ sin 0+. 

l Ml COS 051^ 



+ 4j^|r^+_| sin(Amr + 0+_ - 05W/)e-^^"/2| + ^3^^^ 

(7.6) 

This distribution is less sensitive to a; its leading a effect is of order a ■ \f]^ p. 

However, the measurement of this distribution does allow one to put independent 
constraints on /3 and 7. 

The corresponding distributions for the three pion decay channels can be ob- 
tained by tracing the density matrix P with the operator O^-j^ (8> Ostt as in (6.12). 
Cstt is expressed in terms of the three pion decay amplitudes as 

0„ = |A-3,|^ I f I , (7.7) 
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where 



K2 



(3vr I Ki) 
(Stt I K2)' 



(7.8) 



In this section we ignore CP and CPT violation in the decay amphtudes, so 
we set Y^^ to zero; then 



P(37r,ri;37r,r2) 



X 



,—^STi—TlT2 _j_ g— r^Tl— rsT2 



1 



- 2|%^|2 cos(Am(Ti - r2))e-(^+"-'>')(^i+^^) 

+ ^^\mn\ sin(AmTi - 03^ + 05w)e-(^+a-7)ng-r.r. + (1^2) 



(7.9) 



+ 2 



'7 ,0^1- I '^Stt 



Ar \d\ COS (f)sw. 

where IjysTrle*'^^- = and i?^ = -7/Ar + \fj3n\'^ - 4:{P/Ar)Im[fj3nd/d*] . 
This distribution reduces, on the hne ti = T2 = t, to 



-2rr 



Q(37r;37r; r) oc e 

I COS 0^1^ 



7 



Ar 



+4j^|7737r| sin(AmT - ^stt + 05w)e^^^/^| + 0^{a, P, 7, cm, A). 



(7.10) 



This distribution has a sensitivity to 7 which survives at large time (r^-r 3> 1) 



as 



lim Q(37r; 37r; r) oc -J^e-^^^^ 
rsr»i V ' ' ^ Ar 



(7.11) 



Its origin can be traced to the pl ® PL propagating mode of the density matrix 
(6.7). This behavior contrasts with the behavior of the equal-time distribution 
2(77+77"; tt+tt"; r); the latter vanishes, at large time, as ~ e"^^"^. 
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8. Tests of quantum mechanics at a factory: general final states 



In the previous section, we showed how to constrain the EHNS parameters 
independently of one another by considering the most straightforward experiments 
on decay, analyzed with the simplest theory. In this section, we will generalize 
our analysis to include direct CPT violation and the effects of e'. in the K^^ decay 
matrix elements, as we have discussed for single kaon experiments in Section 5. We 
will first re-examine the determination of the EHNS parameters. We will show that 
these complications have virtually no effect on the method, or even the formulae, 
given in the previous section for the determination of a, /3, and 7. Then we will 
consider the reciprocal problem of the effect of quantum mechanics violation on 
the experimental determination of the kaon decay matrix elements. We will show 
that the measurements of these decay matrix elements can be affected if (5 and 7 
are nonzero. The measurements most sensitive to the modifications of kaon decay 
within quantum mechanics are asymmetries of the integrated time distributions 
(6. 13)!*' We will present formulae which show how these asymmetries are shifted by 
quantum mechanics violation. 

The modification of the formulae for the /; r) in the presence of the effects 
discussed in Section 5 is quite minor. In the formulae for leptonic double decay 
distributions, |ap is changed to |a + 6p for decays to and to |a — 6p for decays 
to ^~ . This has no effect on the functional form of the double time distribution 
except for the simple replacement 



In particular, the formula (7.4) is still valid in this more general context. Similarly, 
the inclusion of more general effects in the kaon decay vertices changes the relative 
normalization of tt+tt" and tt^tt^ decay rates. However, the formula (7.6) remains 
valid with the replacement 



h 

CM + - 

a 




r]+-\e' 



,i4>+ 




(8.2) 
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as in eq. (5.6). This is equally true for the case of the three pion decay, for which 
formula (7.10) remains valid with the substitution* 

Imnle''^'" ^e+ + Ys^ . (8.3) 

Thus, there is no difficulty in constraining CPT violation from outside quantum 
mechanics in (f) factory experiments even in this more general context. For refer- 
ence, the complete expressions for the double time distributions to pion or leptonic 
final states are given in Appendix B. 

It is also interesting to ask whether the converse of this statement is true, 
whether quantum mechanics violation can interfere with the measurement of e' 
and other more standard asymmetries which should appear at a factory. In the 
following, we briefly discuss how the determination of certain of these quantities 
might be affected. 

One would not immediately expect the measured value of e' / e to be significantly 
affected by quantum mechanics violation, since e' is a property of the difference 
between the decay rates into 7r''"7r~ and while the parameters P and 7 affect 
the time evolution of the kaon system prior its decay. Prom this argument, one 
would expect e'/e to be at most corrected by a factor l+0{P/de , 7/|(i||ep). Using 
the bounds on f3 and 7 derived in Section 4, we estimate the terms in parenthesis to 
be at most corrections of the order of 5% and 25% respectively. These corrections 
are mild and can be reduced further if no evidence for nonzero (3 and 7 is found at 
a 4> factory. 

However, if (3 and 7 prove to be nonzero, they will have two important effects 
in the measurement of e'/e. First, any time-dependent method of determining 
e'/e would be complicated by the new time dependences introduced by quantum- 
mechanics violating terms in the density matrix evolution. As we have already 

★ We also assume here that we can neglect a term proportional to ImYa^ ; we expect this term 
to be no bigger than ImY^ Ime'. 
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pointed out, this effect can be minimized by considering the integrated distributions 
(6.13) at fixed time interval Ar = ri — T2. Second, because f3/d and e are close to 
being orthogonal in the complex plane, effects of nonzero /3 can mix the real and 
imaginary parts of e'/e. 

To demonstrate these points quantitatively, we consider the determination'^"' of 
e' / e from the measurement of the quantity 

A (^+^- ■ ^0^0- A.^ = ^(^+^"; At) - Pjn^A vr+TT-; Ar) 

^ ' ' ' :P(7r+7r-; ttOttO; At) + PiTrOn^; n+n-, At) ' ^ ' 

This asymmetry can be straightforwardly computed using the formulae of Ap- 
pendix B. We find 

A(Tr+Tr-; 7r°7r°; Ar) = 3Re - x - 3 Im - x — , (8.5) 

The coefficients //r^i and V are functions of Ar, P/d, 7, and Ifj^ | which are given 

in Appendix C. 

In the context of pure quantum mechanics, the quantities AfRj and V have a 
simple functional form, and the quantities Re e'/e, Ime'/e can be extracted from 
Ar) by a two-parameter fit. In presence of quantum mechanics 
violation, this is no longer true. For example, in the limit F^Ar ^ 1, 

>l(^)(7r+7r-;7rV; Ar) ^ 3Ree7e . (8.6) 

In the same limit, (8.5) becomes 
^(7r+7r-;7rV; Ar) ^ 

l + (2/3/|rf||r?+_|)sin(05W^-0+_) 



3Ree7e 
-3Im e'/e 



1 + (7/AF|^+_|2) + {2p/\d\\f}+-\)ism{2<f>sw - 0+-)/cos05iv) 
{2(3/\d\\f]+-\) cos{(j)sw - 0+- 



(8.7) 

To understand the role of the near-orthogonality oi (3/d and fj^ , we may approx- 
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imate 0-| ~ 4>sw ~ 45°; then (8.7) simplifies to 

^(7r+7r-;7rV; At) ^ 



3Re e'/e 



7 



V2\d\\fj+. 



\V+- 



3 Im e / e 



\V+- 



(8.8) 



This makes clear that, if one opens the possibility of quantum mechanics violation 
in the neutral kaon system at the level /3 ~ m^/mpi, one must constrain or 
determine /3 in order to measure Ree'. 

Similar effects are present in asymmetries involving the semileptonic decay 
distributions P{i^;£^; At). For simplicity, we again consider only the large-time 
limit r^Ar 3> 1. Three particularly informative asymmetries are 



9+ 



At) - V{t 



Ar) + V{t 
Ar) - V{t 



Ar) + V{t 
Ar) - V{i- 



At) 



At) 
At) 



P 



At) 
;Ar) 



4Re (sM + b/a) + 8t-j7 sin 20^^^ cos (l)sw 



\d\ 



V{£+;£+;At) + V{£-]£-]At) 



-4Re A + 4— sin (t)sw [l - 4 cos^ (t>sw\ ■ 



(8.9) 

The first of the above formulae yields a direct determination of 5^.1 even in the 
presence of quantum mechanics violation. However, the other two formulae arc 
more complicated. If quantum mechanics is assumed to be valid, these two limits 
give simple constraints on the CPT violating parameters A and Re(6/a). However, 
in the more general context of quantum mechanics violation, these parameters are 
constrained only to the extent that P is known from the experiments described in 
Section 7. 

It is remarkable how sensitively the parameters of quantum mechanics violation 
affect the various observable quantities of the K^-K^ system as observed at the 
0. Using the strategies we have discussed, it is likely that all three of the EHNS 
parameters of quantum mechanics violation can be bounded, or measured, at a 
level well below (l/mpi). Perhaps there are still more surprises waiting for us in 
neutral kaon physics. 
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APPENDIX A: Eigenmodes of density matrix evolution 



In this appendix, we give more exact formulae for the eigenmodes of the density 
matrix evolution and for the components of the density matrix for neutral kaon 
pairs resulting from decay. The expressions below are complete through second 
order in CP- violating parameters. 

First, we present the density matrix eigenmodes. These are given as column 
vectors 

P2 
PI 
\Pl/ 



(A.l) 



where these elements are defined by eq. (3.13). The expressions (3.22) should be 
replaced by: 



PS 



( 1 \ 

-7/Ar + 2Im[e+e^*(i]/Ar 

e^* - (cie+7/d*Ar) + i(e^*7/d*) - (2Q;Ime^]/d*) 
V - (rf*e;^*7/dAr) - i[Cg-ild) - (2alme^]/(i) / 



(A.2) 



PL 



7/Ar - 2Im[e+e^*d*]/Ar 
1 

4 + (rf*es*7/c^Ar) + 2(6+7/ri) + (2aIme+]/ri) 
Ve+* + (de^7/o?*Ar) - i(e+*7/o?*) + (2aIme+]/d*) / 



(A.3) 
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PI 



( e^^*{l + i-i/d* - ia/d*) - ie~^j/d* + i{de]^a/2d*Am) \ 
ej(l + i-f/d - ia/d) - ie2*-f/d + i(d*epa/2dAm) 

1 



V 



(A.4) 



-i{a/2Am) + {tgt-^*d + e+*e+*rf*)/2Am j 



/ e;^(l - i-^/d + ia/d) + iep-f/d* - i{d*t^*a/2dAm) \ 
e+*(l - i-f/d* + ia/d*) + ie;^7/d* - i{de^a/2d*Am) 



Pi 



i{a/2Am) + {e^*eld* + e+eJ*(i)/2A 



m 



(A.5) 



Note that the expressions for pLi and ps2 given in (3.22) agree with the expressions 
above. 



Next, we present a more precise form for the density matrix P by quoting the 
matrix elements of Aij, defined by eq. (6.11), to second order in small parameters. 
We find: 
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2' V o . ^ |rf|4 2 (Ar)2 

A ^ ^4 T 

Asi ^ Ais ^ (A^jY ^ {AjsT 

2 - ImcL - z - (ei - e^) + z ■ 



d ^ " 2Amd* AT d* 

Ali = ^/L = {AaT = (^7l)* 



d* d* ' 2Amd AFd 

^// = An = (1 + 2ReeLe5) + /?Re 

3 2 ci^ + (rf*)^ _ 3 
^2^ 8 (A^ 

yl^^ = (^^)* = _ 2 ^ - (ei - e^) 

^ 2Am MP Am ^ 



(A.6) 



APPENDIX B: Double time distributions 



In this appendix, we give the complete formulae for the double time distribu- 
tions V{fi,Ti; /2,T2) for (p decay to the various final states discussed in the text. 
The expressions below are complete through first order in CP- violating parameters 
unless it is specified otherwise. 
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7'(£±,ri;£±,r2) = ^ 
I [(1 ± 4Re (eM + ^)] [e"^^^^"^^^^ + "^^^'^ '''^ - 2 cos(Am(ri - r2))e-(^+"-^)(^^+^^)] 
± 4^ sin(Amri - (j)sw)e~^^^'^~^^^'e-^'^^ + (1^2) 

± 4^sin(AmTi + (t)sw)e~^^^'^~^^^'e-^'^^^ + (1^2) 

\a\ 

+ 2-^ sin Amfri + ra) e-(^^+"-T)(^i+^2) + 2^ rg-r^n+ra) _ g-rs(n+r2)j \ 
Am Ar '- J 

(B.l) 



X 



7'(£+,ri;r,r2) = ^ 

X I (1 + 4Re (A - /?/ci))e-^^^i-^^""^ + (1 - 4Re (A - {3 / d))e-^ '^^''^ 
+ 2 cos(Am(ri — r2 ) - 4Im(A + (3/d)) e-(r+«-7)(ri+r2) 
+ 4^sin(AmTi - 05^^)6" ^^+""^^^^6-^^^^ -(1^2) 

+ 4^sin(Amri + 4)swy~^^^'^~^^^'e-^'^^^ -(1^2) 
\a\ 

_ 2-^ sin Amfri + ra) e-(^+«-7)(Ti+r2) ^ re-r4n+^2) _ g-rs(ri+r2)i \ 
Am ^ ' Ar^ -'J 

(B.2) 



x|i?L[e-^«^i-^^^2 + e-r^^i-^s^^] - 2|?7+_|2cos(Am(Ti - T2)) e-(^+«-7)(n+r2) 

+ 44|?7+_| sin(AmTi + - 05Ty)e-^^+"-'^)^^e-r^^^ + (1^2) 
\d\ 

-2(^ + 4Aini[^+_-r+_]^e-r«(-^+-^)| + ^^(a, ^, e^^^, , 

(B.3) 
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using ?7_| 6**^+" = + and 

Rl = 7/Ar + |?7+_|2 + 4(/3/Ar)Im[?7+_d/d* - . (B.4) 



^/ + - \ l-'^H—Pl-^OOp 
V^-K^-K ,Ti;7r TT ,T2) = ^ 

- 2|77+_||77oo| cos(Am(ri - ra) + 0+- - 0oo) e-(r+«-T')(^i+^^) 
+ 4^|r/+-| sin(AmTi - (t>sw + 0+-)e-(^+«-7)^ie-rs^2 

+ 4y||r7oo| sin(AmT2 - ct>sw + ^ooje-^^+'^-^^'^e-r^"^ 

-2(^^ + 2^Im[r/+_-F+_ + 77oo->oo])e-r^-(^^+^^)| + 0\a, (3,^,es,L.Y^ 

(B.5) 



P(37r,Ti;37r,r2) = — ^ 

- 2 1 77377 P cos ( Am (n - r2))e-(^+'^-^)(^i+^^) 

+ 4^|r/3.| sin(Amri - 03, + 0sw^)e-(r+«-T)-ie-r-^^ + (1^2) 

+ 2(^^ + 4^Im[7737r-l3,.])e-r-(^^+^^)| + 0\a,(3,^,es^L,Y:,,). 

(B.6) 

where |r/37r|e*''^=^'' = + -^Stt 

i?^ = -7/Ar + |?73^|2 - 4(/3/Ar)Im[7?3^ci/d* - Y^^] . (B.7) 

45 



P(7r+7r-,ri;£±,r2) 



|X+_|2|a|2 



-FsTl— rz,T2 



^ 2|?7+_| cos(Am(Ti - r2) + 0+_)e-(^+"-T)(^i+^^) 
± 4^ sin(AmT2 - (f)sw)e~^^+''~^^^'e-'^'^' 



3-rs(Tl+T2) 



(B.8) 



APPENDIX C: Formulae for the measurement of e'/e 

In this Appendix, wc present formulae for extracting e'/e from the integrated 
distributions at fixed time interval At = t\ — T2 of the asymmetric decay into 
charged and neutral pion final states. 



^(7r+7r-;7rV;AT) = 



P(7r+7r-; ttV; At) - P(7rV; tt+tt"; At) 
P(7r+7r-; ttOttO; At) + V{t^^t^^\ tt+tt"; At) ' 



(C.l) 



>l(7r+7r-;7rV; At) = SRee'/e x ^ - Slme'/e x 



with 



(C.2) 



1 + 2 



■sin(05;^ - 0+-) 



1 + ^ i^ii^^ | (sin(0giy - 0+-) - 1^1 sin((/)5VF + 0« - 0+-)) 



+e 



-FAt 



/3 



1^1 sin(AmAT + 0+_ — 4>sw — 4>z) 



(C.3) 
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Ml = e-r-^- 



/3 


\d\ 


l^+-l 






\d\ 





cos{(i)sw - (t>+-) 



■)) 



+e 



-rAr 



2 sin Am At - 2- 



l^+-l 



1 2; I cos(AmAT + 0-1 — ~ 



(C.4) 



-rsAr 



+e 



— e 



1 + 
1 - 



7 



+ 2, 



Ar|77+_|2 
Ar|fj+_|2rs 



+ 2 



_^ sm(20gH^ - 0+_ 

|r;+_| COS05VK 

/? fsui{2(f)sw - </>+-) 



COS (l)sw 



- 2\z\ s\n{(t)sw + (pz- 0+-)) 



-rAr 



2 COS Am At — 4- 



/3 



d\\rj+. 



\z\ sin(AmAT + 0-| — (psw 



W t have defined 



2r 



+ r + iAm 



(C.5) 



(C.6) 
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FIGURE CAPTIONS 

1) Constraints on the EHNS parameters P and 7 from the comparison of deter- 
minations of the e parameter from different observables of the K^-Ks system: 
(a) the systematics of expected discrepancies; (b) the current experimental 
situation. 
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